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Abstract 

The present experimental status concerning the composition of cos- 
mic rays in the PeV region is reviewed. The results are compared to 
predictions of contemporary models for the acceleration and propaga- 
tion of galactic cosmic rays. 



1 Introduction 

The solar system is permanently exposed to a vast flux of highly energetic 
and fully ionized atomic nuclei, the cosmic rays. Their energies extend from 
the GeV range to at least 10^" eV. Over a wide range the energy spectrum 
follows a power law dN/dE oc E'^ . The spectral index changes around 4 PeV 
from 7 ~ —2.7 to 7 ~ —3.1. This transition phenomenon has been reported 
for the first time in 1958 and is commonly referred to as the knee. A 
cutout of the energy spectrum in the knee region is depicted in Fig.^ Shown 
is the all-particle energy spectrum, which exhibits the knee, as well as results 
from direct measurements above the atmosphere for primary protons and 
iron nuclei. 

The origin of the knee is one of the central questions of high-energy as- 
troparticle physics, closely related to the mechanisms of acceleration and 
propagation of high-energy cosmic rays. Answers are expected from the 
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Figure 1: Energy spectrum of 
cosmic rays for all particles (dark 
grey), protons (black), and iron 
lo' lo' lo' 10= 10' 10' 10= 10' nuclei (light grey), for references 

Energy E„ [GeV] gge [Jj. 

measurements of the energy spectra for individual elements (or at least ele- 
mental groups) above 1 PeV. Since this is an extreme experimental challenge, 
one often measures the average atomic mass instead. 

In this article, the present experimental status is reviewed and the results 
are set in context to predictions of contemporary models for the acceleration 
and propagation of cosmic rays. 

Ideally, one would like to continue the energy spectra for individual el- 
ements above several 10^^ eV with direct measurements above the atmo- 
sphere. At present, several groups measure in this region, applying different 
experimental techniques, e.g. a calorimeter (ATIC jH]) or a transition ra- 
diation detector (TRACER 4 ). The TRACER experiment, investigating 
nuclei from oxygen to iron, had a 14 day flight from Mc Murdo, Antarctica 
in December 2003. Recent accelerator tests have shown that transition ra- 
diation detectors can be utilized to measure energy spectra of cosmic rays 
in a space borne experiment up to energies of about 1 PeV/n 0. 

Presently, at energies above 1 PeV one relies on indirect measurements. 
In these experiments the secondary products, generated by interactions of 
cosmic-ray particles with nuclei in the atmosphere are investigated. Two 
basic approaches can be distinguished: Measuring the debris of the par- 
ticle cascade at ground level by registering electrons, muons, or hadrons. 
Or measuring the longitudinal shower development in the atmosphere by 
exploration of the Cerenkov or Fluorescence light generated predominantly 
by the shower electrons. An astrophysical interpretation of air shower data 
requires detailed knowledge of the interaction processes in the atmosphere. 
One of the tasks of air shower experiments is to improve the understanding 
of high-energy interactions above the energies covered by todays accelerator 
experiments and beyond their kinematical bound, see e.g. 




2 Origin of the knee 



The bulk of cosmic rays is assumed to be accelerated in shock fronts of su- 
pernova remnants (SNRs). This goes back to an idea of Baade and Zwicky, 
who have estimated the total power required to generate the observed flux 
of cosmic rays (Tj. It can be shown, that about 3 supernovae per galaxy and 
century are sufficient to release enough kinetic energy in order to deliver 
the required power. A mechanism to accelerate particles by moving mag- 
netic clouds has been introduced by Fermi [8;. The present understanding 
of acceleration in strong shock fronts has been initiated by Blanford and 
Ostriker |2j which could demonstrate that at strong shocks particles are ac- 
celerated efficiently. The finite lifetime of a shock front (~ 10^ a) limits 
the maximum energy attainable to Emax ~ Z ■ (0.1 — 5) PeV for particles 
with charge Z. Various versions of this scenario have been discussed, see 
e.g. ^l^n^lElE]- In the literature also other possibilities, like the 
acceleration of particles in 7-ray bursts are discussed ^3 1^ ■ 

After acceleration, the particles propagate in a diffusive process for about 
20 • 10® a through the Galaxy, being deflected many times by the randomly 
oriented galactic magnetic fields ~ 3 /iG). The propagation is accompa- 
nied by leakage of particles out of the Galaxy. With increasing energy it is 
more and more difficult to magnetically bind the nuclei to the Galaxy. The 
pathlength of traversed material decreases as A oc , with 6 ~ —0.6. Many 
approaches have been undertaken to describe the propagation process, see 
^1 ^1 El EOl 1^ • During the propagation phase, reacceleration of particles 
has been suggested at shock fronts originating from the galactic wind j22j . 

Both, the acceleration and the propagation/leakage processes are ex- 
pected to yield cut-offs for the fluxes of nuclei at energies proportional to 
their nuclear charge E^ oc Z. 

As further possible causes for the knee, interactions of cosmic-ray nuclei 
with background photons or neutrinos in the Galaxy [23L 1241 OS] or new 
types of interactions in the atmosphere [SHI have been discussed. For a 
more comprehensive review of models the reader is referred to j27j . 

Electromagnetic emission of SNR has been detected in a wide energy 
range from radio wave lengths to the x-ray regime. The observations can be 
interpreted as synchrotron emission from electrons, which are accelerated 
in these regions j29j . Recently, the HEGRA experiments has detected an 
excess of high-energy 7-rays from the supernova remnant Gas A, see Fig.|21 
This is interpreted as evidence for hadron acceleration in the SNR. The 
hadrons interact with protons of the interstellar medium, producing vr'^s, 
which decay into high-energy photons, supposedly detected by the HEGRA 




Figure 2: Left: The dots indicate the measured 7-ray signal as function of angular 
distance to Cassiopeia A as measured by the HEGRA experiment. The shaded 
area is a background estimate |28| . Right: Integral 7-ray fluxes due to inverse 
Compton scattering, non-thermal bremsstrahlung and 7r°-decay according to the 
model by Berezhko et al. compared with measurements of the HEGRA and EGRET 
experiments j^lj. The dotted line is obtained for an assumed cut-off at 4 TeV. 

experiment j29j . The flux is compatible with a model of electron and hadron 
acceleration in shock fronts, see Fig.|2l 

Recently, an excess of charged particles from the direction of a SNR 
(Monogem ring, d ~ 300 pc) has been reported [30J . However, such a signal 
could not be confirmed by the KASCADE experiment |3]j|. The latter has 
performed a detailed search for point sources, covering the whole visible sky 
at energies Eq > 0.3 PeV. Special attention has been given to the region of 
the galactic plane, as well as to the vicinity of known SNRs and TeV-7-ray 
sources. No significant excess could be found. The gyromagnetic radius^ 
of particles with an energy around 1 PeV in the galactic magnetic field is 
in the order of 1 pc. Hence, it is not expected to find any point sources at 
these energies. 

To characterize the large scale anisotropy of cosmic rays, the Rayleigh 
amplitude 
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\ 




for n events which are detected at a right ascension angle is introduced. A 
compilation of measured amplitudes as function of energy is given in Fig.|51 
j32j . An increase of A as function of energy can be recognized. This trend 

^ro — p/(ZeB) for a particle with momentum p and charge Ze in a magnetic field B. 
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Figure 3: Rayleigh amplitudes as 
function of energy for various exper- 
iments, for references see j^. Addi- 
tionally, model predictions for a dif- 
fusion model are shown. The lines 
indicate the expected anisotropy for 
primary protons (H) and iron nuclei 
(Fe) as well as for all particles (tot) 
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is compatible with expectations taking into account diffusive propagation of 
cosmic rays in the Galaxy as indicated by the lines. 



3 Energy spectrum of cosmic rays 

The all-particle energy spectrum has been derived by many experiments as 
shown in Fig.|ll In this representation the flux of the air shower experi- 
ments has been normalized to the extrapolated flux of direct measurements 
at 1 PeV by introducing a slight adjustment (±10%) of the energy scales 
of the individual experiments '7 . It is quite interesting to realize that the 
absolute energy calibration of these various experiments, using different ob- 
servation techniques and models to describe the shower development in the 
atmosphere agree within ±10%. The normalized all-particle flux changes 
smoothly without any prominent structures. The solid and doted lines in- 
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Figure 5: Results from the KASCADE experiment: Energy spectra for groups 
of elements derived from the data using CORSIKA with the hadronic interaction 
models QGSJET (left) and SIBYLL (right) ||34j. In addition, the proton spectrum 
as derived from the measurements of single hadrons is shown • The lines indicate 
extrapolations of direct measurements [2]- 

dicate the total sum of galactic cosmic rays according to a parameterization 
of data _2 . 

The KASCADE group performed systematic studies to evaluate the in- 
fluence of different hadronic interaction models used in the simulations to in- 
terpret the data on the resulting spectra for elemental groups [HH ■ Two sets 
of spectra, derived from the observation of the electromagnetic and muonic 
air shower components, applying the Gold algorithm and using CORSIKA 
|37j with the hadronic interaction models QGSJET and SIBYLL are com- 
piled in Figure |SJ 

As can be seen in the figure, the flux for elemental groups depends on the 
model used. The KASCADE group emphasizes that at present, there are 
systematic differences between measured and simulated observables, which 
result in the ambiguities of the spectra for elemental groups. These conclu- 
sions apply in a similar way also to other experiments. A correct deconvo- 
lution of energy spectra for elemental groups requires a precise knowledge 
of the hadronic interactions in the atmosphere. The interaction models 
presently used do not describe the measurements with a sufficiently high 
precision required for this task. 

The figure also shows the spectrum of primary protons, which has been 
derived from the flux of unaccompanied hadrons jSS]. The spectrum is com- 
patible with the proton flux as obtained from the unfolding procedure when 
using the QGSJET model. For comparison, also the results of direct mea- 
surements at lower energies at the top of the atmosphere [2] are presented 
in the figure. 
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Figure 6: Energy spectra for elemen- 
tal groups. Top left: protons, top 
right: helium, and bottom: iron nu- 
clei. Open symbols give results of di- 
rect measurements, for references see 
Filled symbols represent results from air 
shower measurements: KASCADE elec- 
trons/muons interpreted with two inter- 
action models (preliminary), KAS- 
CADE single hadrons and EAS-Top 
The data are compared to predic- 
tions of calculations by Kalmykov et al. 
PS] (•••), Sveshnikova et al. [El (---), 
and the Poly-Gonato model (2] ( — ). 



In order to give an impression of the present status, the energy spectra 
for three elemental groups (protons, helium, and iron nuclei) are compiled 
in Fig. ini over a wide energy range. Shown are results from direct measure- 
ments above the atmosphere as well as the KASCADE results described 
above. The EAS-TOP experiment published two sets of spectra with dif- 
ferent assumptions about the contribution of protons and helium nuclei, for 
details see |H^. The resulting fluxes are indicated by two squares per pri- 
mary energy. To guide the eye, the solid lines indicate power law spectra 
with a cut-off at Z ■ 4.5 PeV. For the iron spectrum at low energies the 
influence of modulation due to the magnetic fields of the heliosphere can be 
recognized. 

The dashed lines represent calculations of energy spectra for nuclei accel- 
erated in supernova remnants by Sveshnikova et al. 13 . It is assumed that 
the particles are accelerated in a variety of supernovae populations, each 
having an individual maximum energy to be attained during acceleration, 
which results in the bumpy structure of the obtained spectra. The dotted 



lines reflect calculations of the diffusive propagation of particles through the 
Galaxy by Kalmykov et al. [^Hl- The leakage of particles yields a rigidity 
dependent cut-off. Comparison with the data may suggest a qualitative 
understanding of the energy spectra. However, for a precise quantitative 
understanding, detailed investigations of the systematic errors of the mea- 
surements are necessary and the description of the interaction processes in 
the atmosphere needs to be improved. 

4 Mass composition of cosmic rays 

The position of the maximum of an electromagnetic or nuclear cascade in 
matter depends on the incident particles energy as X^ax oc In En, where En 
is the energy per nucleon En = Eq/A. Hence, Xmax depends on In^. Also, 
many other observables in air showers like the number of electrons, muons, 
or hadrons observed at ground level depend roughly on In A. To characterize 
the cosmic-ray mass composition one uses commonly the mean logarithmic 
mass (In A) , defined as (In A) = X) ^« 1^ ) where is the relative fraction 
of nuclei with atomic mass number Ai. 

The mean logarithmic mass from experiments measuring electrons, 
muons, and hadrons at ground level using mostly CORSIKA/QGSJET to 
interpret the data are compiled in Fig.[7l;. One recognizes an increase of 
(In A) as function of energy, the line indicates an increase as expected for a 
rigidity dependent cut-off for individual elements HJ. 

A second group of experiments measures the average depth of the shower 
maximum Xmax in the atmosphere by registering Cerenkov photons or fiuo- 
rescence light. Interpreting this data with CORSIKA/QGSJET 01 leads to 
the results shown in Fig.[7K. An increase of (In A) as described before is not 
obtained. 

The description of hadronic interactions within the models is based on 
accelerator measurements. However, extrapolating the experimental uncer- 
tainties to higher energies may cause significant differences in the interpre- 
tation of air shower data. For example, the cross sections for proton - anti 
proton collisions at the Tevatron (y/s = 1.8 TeV) exhibit an error of about 
10%. To investigate the influence of such uncertainties on the development 
of air showers, in the model QGSJET 01 the cross section for p-p collisions 
has been reduced to the lower error boundary of the Tevatron data. This 
results in a deeper penetration of showers into the atmosphere. Accordingly, 
{In A) derived with such an interaction model, using the same experimen- 
tal Xjnax-values as before, yields larger values at high energies, as shown 
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Figure 7: Mean logarithmic mass as function of energy as obtained by various 
experiments. Top: (In^) derived from observations of the average depth of the 
shower maximum Xmax interpreted with CORSIKA/QGSJETOl (a) and with a 
modified version (b) j39J . Bottom: Experiments measuring electrons, muons, and 
hadrons at ground level. For details and references see pi 139). 
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Figure 8: Mean logarithmic mass as function of energy as obtained by various 
experiments (shaded area) compared with different models (lines), a) Acceleration 
in SNR (Berezhko et al. ^U] i Stanev et al. , Kobayakawa et al. ^21 1 Sveshnikova 
et al. 5S|); b) Single source model (Eriykin & Wolfendale ^1]), acceleration in 
7-ray bursts (Plaga Jl], Wick et al. ^El), reacceleration in the galactic wind (Volk 
et al. 1221); c) Diffusion in Galaxy (Kalmykov et al. Ogio et al. Roulet 
et al. ^ni); d) Propagation in the Galaxy (Lagutin et al. |2I], Swordy ^3) and 
interaction with background photons (Tkaczyk |23j) and neutrinos (Dova et al. 
121). For details see 127]. 

in Fig.[7|3 [Sni. These values follow better the trend of the observations of 
particles at the ground, shown in Fig.[2h- 

Combining the experimental values from Fig.[7|3 and c, the result is dis- 
played as grey band in Fig.|SJ The experimental data are going to be com- 
pared to various models, which are discussed in the literature as possible 
origin of the knee j27j . 

Several approaches to model the shock acceleration in SNRs (Berezhko 
et al. ) Kobayakawa et al. ^21 ) Stanev et al. , Sveshnikova et al. ^S] ) 
are summarized in Fig.|HK. The models differ in assumptions of properties 
of the SNRs like magnetic field strength, available energy etc. This yields 
differences in (In A), as can be inferred from the figure. 



Predictions of the single source model (Erlykin & Wolfendale ^1] ) , reac- 
celeration in the galactic wind (Volk et al. [22]) and acceleration in 7-ray 
bursts (Plaga Wick et al. ^Hl) are shown in Fig.|Hl3. The latter differ 
in their interpretation of the origin for the knee, Plaga attributes it to the 
maximum energy reached during the acceleration process, while Wick et al. 
propose leakage from the Galaxy as the cause. 

The propagation of cosmic rays as described in diffusion models 
(Kalmykov et al. |HH1, Ogio et al. [lH|, Roulet et al. yields (In^)- 

values presented in Fig.|St. They are all based on the approach by Ptuskin 
et al. ^Zj, but take into account different assumptions on details of the 
propagation process, like the structure of galactic magnetic fields etc. 

The last panel (Fig.|Hll) summarizes predictions of models taking into 
account cosmic-ray propagation in the Galaxy (Lagutin Swordy 20 ) as 
well as interactions with background photons (Tkaczyk ) and neutrinos 
(Dova et al. [21). 

As can be inferred from Fig.[Hl the situation is not yet conclusive, but a 
trend towards a standard picture can be recognized. Some of the proposed 
explanations can already be excluded. Interactions with background parti- 
cles in the Galaxy would produce a big amount of secondary protons, which 
results in a light mass composition at high energies, not confirmed by the 
experiments. Furthermore, a massive neutrino, proposed in |^ is excluded 
by measurements of the WMAP and 2dFGRS experiments [in|. The ap- 
proach described in (T5l does not describe the trend of the data. The data 
are consistent with acceleration in SNRs (Fig.[HK) and diffusive propagation 
(Fig.[Hl;). Also reacceleration in the galactic wind (Volk et al.) and accel- 
eration in 7-ray bursts in combination with diffusive propagation (Wick et 
al.) follow the trend indicated by the data. 

5 Conclusions 

Comparing the present results to the status before the start of the KAS- 
CADE experiment (about one decade ago), our knowledge about high- 
energy cosmic rays has significantly improved. The experiment has shown 
that the knee is caused by the subsequent cut-offs of individual elements, 
starting with protons and helium nuclei and that the mean logarithmic mass 
increases as function of energy. 

Summarizing the large number of experimental observations, there are 
indications for a standard picture of galactic cosmic rays. At least a large 
fraction of them seems to be accelerated in supernova remnants up to ener- 



gies of Z • (0.1 — 5) PeV. Higher energies may be reached by reacceleration 
in the galactic wind or by acceleration in additional sources, such as 7-ray 
bursts. The particles propagate in a diffusive process through the Galaxy. 
With rising energy the pathlength decreases and particles escape easier from 
the Galaxy. This brings about the knee in the energy spectrum. 

When the understanding of the hadronic interactions in the atmosphere 
improves, the measurements can be interpreted with higher reliability. This 
will put more restrictions on the models to describe acceleration and prop- 
agation of cosmic rays. 
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